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“ J  A  formalism  is  developed  for  calculating  the  complex  refractive 

■  index  of  a  monatomic  gas  as  a  function  of  frequency  in  the. region  of  an 

i"~~^  absorption.  The  treatment  involves  a  direct  study  of  the  properties  of 

U_l 

CJD  the  dressed  photon  state  in  the  gas,  and  includes  effects  due  to  trans- 

e'--C  lational  motion  of  the  atoms  as  well  as  the  dipolar  "resonance"  interac- 
tions  between  them.  It  is  necessary  to  assume  that  the  electromagnetic 


CO 
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properties  of  the  gas  are  describable  in  some  detail  by  a  linear  frequency 

and  wave  number -dependent  susceptibility.  A  set  of  coupled  nonlinear 

integral  equations  are  derived  which  together  determine  the  susceptibility 

function  and  hence  the  observable  refractive  index.  The  "static"  limit  of 

large  atomic  mass  is  considered  in  some  detail,  and  a  first  correction  to 

it  is  also  obtained.  The  results  are  compared  with  measurements  by 

Tomiser  on  the  linewidths  of  the  sodium  D-line  at  various  temperatures 

and  pressures.  We  obtain  qualitative  agreement,  while  previous  theoret«- 
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ical  linewidths  were  too  small  by  a  factor  of  order  10  . 
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I.  INTRODUCTION 

The  problem  of  calculating  absorption  line  shapes  in  monatomic 
gases,  especially  the  contribution  of  the  dipolar  "resonance"  interactions 
between  different  absorbers,  has  been  attacked  in  the  past  in  several 
different  ways,  all  of  a  more  or  less  ad  hoc  nature?  The  various  methods 
("collision",  "statistical")  are  in  agreement  among  themselves  at  least 
as  to  the  qualitative  features  of  the  line  shape,  and  as  to  the  order  of 
magnitude  of  the  linewidth.  However,  it  appears  to  have  gone  largely 

unnoticed  that  they  are  in  strong  disagreement  with  the  latest  experimental 

2,3  1 

results.  In  particular,  the  calculated  "resonance"  width  is  too  small 

3 

by  a  factor  of  the  order  of  10  to  fit  the  experimental  data,  and  the  Doppler 
width  is  also  much  too  small  to  account  for  the  observed  effects.  (Moser 

3 

and  Schultz  give  a  detailed  comparison  of  the  theoretical  results  with 

their  own  and  other  experiments. ) 

4  5 

One  of  us  ’  has  developed  a  theory  of  the  line  shape  which  begins 
from  the  more  fundamental  point  of  view  of  actually  studying  the  properties 
of  the  stationary  states  of  the  quantum  mechanical  system  composed  of 
radiation  field  and  matter.  This  theory,  however,  considered  each  absorber 
to  occupy  a  fixed  position  in  space,  so  that  effects  due  to  the  motion  of  the 
atoms  were  left  out.  The  original  formulation  depended  on  the  use  of 


R.  G.  Breene,  Jr.,  Revs.  Modern  Phys.  29,  94  (1957) 

^  J.  Tomiser,  Acta  Phys.  Austrlaca  8,  198  (1953);  Q,  276  (1954);  8,  323 
(1954);  J,  18  (1954). 

3  H.  Moser  and  H.  Schultz,  Ann.  Physik_7,  243  (1959). 

4  C.  A.  Mead,  Phys.  Rev.  120,  854  (I960). 

5  C.  A.  Mead,  Phys.  Rev.  120,  860  (I960). 
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"dsmping  operators"  ,  but  more  recently  the  results  have  been  rederlved 
by  means  of  a  simpler  formalism,  and  at  the  same  time  brought  into  rough 
agreement  with  experiment.  ^ 

The  purpose  of  the  present  article  is  to  remove  the  defect  of  neglecting 
the  translational  motion  of  the  absorbers,  and  also  to  discuss  in  somewhat 
more  detail  than  before  the  comparison  of  theory  with  experiment.  Our 
method  will  be  a  generalization  of  that  of  reference  7;  that  is,  we  will 
study  the  properties  of  the  stationary  states  of  the  system  but  will  also 
assume  that  the  system  is  deacribable  by  means  of  a  refractive  index. 

This  assumption  will  be  generalized  to  take  account  of  the  possibility  of 

g 

"spatial  dispersion".  The  plan  of  the  article  is  as  follows:  In  section 

II,  the  notation  is  introduced  and  some  preliminary  matters  discussed, 
including  the  "cutoff  function"  inserted  in  certain  Hamiltonian  matrix 
elements  which  is  supposed  to  represent  roughly  the  effect  of  finite  atomic 
size  and  which  later  prevents  certain  integrals  from  diverging.  In  section 

III,  the  basic  assumption  that  the  medium  is  describable  in  terms  of  a 
refractive  index  is  stated  precisely  as  a  statement  about  the  stationary 
states  of  the  system.  Section  IV  consists  of  the  derivation  of  a  set  of 
integral  equations  which  determine  the  refractive  index.  Section  V  dis¬ 
cusses  the  "static  limit"  of  large  atomic  mass  (which  corresponds  to  the 
situation  in  references 4,  5,  and  7)  and  also  the  lowest  order  correction  for 

finite  mass.  It  is  shown  that  under  the  conditions  of  Tomlser's  experiments 
2 

on  Sodium  the  static  approximation  should  be  adequate.  Section  VI  gives 
a  more  detailed  comparison  with  Tomiser's  results  in  the  static  limit, 

6  C.  A.  Mead,  Phys.  Rev.  112,  1843  (1958). 

7  C.  A.  Mead,  Phys.  Rev.  128,  1753  (1962). 

8  J.  Neufeld,  Phys.  Rev.  123,  1  (1961). 
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including  the  effect  of  the  spin-orbit  splitting.  We  give  a  detailed 
comparison  only  with  Tomiser's  experiments  on  the  Sodium  D-line;  the 
situation  with  the  other  experiments  of  references  2  and  3  is  qualitatively 
similar.  A  brief  discussion  of  the  cutoff  function  used  is  given  in  an 
appendix. 

A  subsequent  paper  by  one  of  us  will  discuss  another  method  for 
obtaining  values  for  certain  integral  properties  of  the  line  shape,  including 
corrections  for  the  failure  of  one  of  the  assumptions  of  Sec.  III. 

II.  NOTATION  AND  HAMILTONIAN 

Our  system  consits  of  N  identical  absorbers,  each  with  mechanical 
mass  m,  distributed  in  a  volume  V,  within  which  the  radiation  field  is 
required  formally  to  obey  periodic  boundary  conditions.  In  the  "initial" 
state  0>,  each  atom  A  has  a  momentum  assigned  to  it,  and  the 
number  of  atoms  with  momentum  in  a  given  region  of  wave-number  space 
is  described  by  a  Boltzmann  distribution  at  some  temperature  T.  For  the 
sake  of  definiteness,  the  ground  electronic  state  of  each  atom  is  considered 
to  be  a  nondegenerate  s -state,  and  the  only  excited  level  (or  at  least  the 
only  one  of  interest  to  us)  is  assumed  to  be  a  triply  degenerate  p-level. 

The  energy  difference  between  the  two  levels  is  |i  vQ.  A  given  state  of 
the  system  of  radiation  field  plus  matter  may  differ  from  the  state  |  0^>  in 
one  or  more  of  three  ways:  (1)  There  may  be  a  photon  present;  this  is 
denoted  by  a  greek  letter,  e.g.  X.  denotes  the  presence  of  a  bare  photon  of 
wave  number/t  .  and  polarization  .  (2)  One  or  more  of  the  atoms  may 

be  excited;  this  is  indicated  by  a  capital  Latin  letter  denoting  which  atom  is 
excited,  with  a  subscript  giving  the  direction  of  "polarization"  (electric 
dipole  moment  matrix  element  between  excited  state  and  ground  state). 

Thus  A^  represents  the  excitation  of  atom  A,  with  dipole  moment  matrix 
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element  in  the  i-direction  (i  =  x,  y,z).  (3)  One  or  more  of  the  atoms  may 

have  acquired  a  momentum  different  from  that  assigned  to  it  in  the  state 
|  0^  .  K,  (A)  indicates  that  atom  A  has  an  excess  momenturn'fi^  ,  so 
that  its  total  momentum  is  )•  The  way  in  which  these  symbols 

are  combined  in  specifying  a  state  is  best  explained  by  means  of  a  few 
examples  of  states  represented  in  our  notation  and  with  their  properties 
also  written  out  in  words: 


I 


’  Photon  present  with  wave  number  A  ^  and  polar¬ 
ization  £  m  ;  atom  A  has  excess  momentum  ft; 
all  other  atoms  as  in  state  |  0  ^  . 


(a1^(A)  Atom  A  excited,  with  "pblarization"  in  i-direction; 

atom  A  has  excess  momentum ;  atom  B  has 
excess  of all  other  atoms  as  in  state  |  0^  ; 
no  photons  present. 


We  are  interested  in  constructing  the  "dressed"  photon  state 
I  X^,  which  is  an  eigenstate  of  the  complete  Hamiltonian  with  energy 
v.  (v  =  c  ^  ),  and  which  in  some  sense  approaches  the  bare  photon 
state  |  X>  in  the  limit  of  zero  interaction.  We  assume  that  |  X  ^  can 
be  built  up  entirely  out  of  two  kinds  of  states: 

(1)  States  with  one  photon  present,  and  no  electronic  excitation. 

(2)  States  with  one  and  only  one  atom  excited,  and  no  photons 
present. 

States  coming  under  categories  (1)  and  (2)  may  still,  of  course,  include 
any  combination  of  excess  translational  momenta  for  the  various  atoms. 
Other  states,  such  as  two-photon  states  and  states  with  several  atoms 
excited,  may  be  included  later  by  means  of  perturbation  theory,  if  desired. 
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The  approximation  of  limiting  ourselves  to  the  states  of  types  (1)  and  (2) 

should  hold  in  the  frequency  region  v.  *=  v  ,  which  is  the  region  we  are 

\  o 

interested  in  here.  It  could  be  made  more  formally  precise  by  means  of 

4  9  10 

the  Arnous -Bleuler  transformation,  ’  ’  but  this  would  make  some  of 
the  subsequent  manipulations  more  complicated  so  it  was  decided  not  to 
use  it  in  this  article. 

9  E.  iirnous  and  K.  Bleuler,  Helv.  Phys.  t.  eta  25,  581  (1952). 

10  W.  Heitler,  The  Quantum  Theory  of  Radiation  (Oxford  University  Press, 
New  York,  1954),  3'rd  ed,  pages  348-353. 


The  refractive  index  can  be  thought  of  as  describing  the  propagation 
of  the  electric  field  matrix  element 

<•!£.<£>  I  *>  *  (5,)^><t|r>- 

This  immediately  leads  to  a  difficulty:  It  would  appear  that  because  of 
momentum  conservation^  I  =  0  unless  =  /l  ^  •  This  would 
mean  that 

<°l  Ea  (R)  I  ^)>  =  £0 oexp(< 

corresponding  to  a  refractive  index  of  unity.  This  problem  does  not  arise 
in  the  situation  treated  in  references  4,  5,  and  7,  since  there  the  momentum 
conservation  is  spoiled  by  holding  the  atoms  fixed.  This  apparently  hope¬ 
less  difficulty  really  is  the  result  of  a  rather  trivial  effect:  Momentum  is 
indeed  conserved  when  an  electromagnetic  wave  propagates  through  a 
medium  with  altered  wave  number,  but  the  excess  momentum  is  taken  up 
either  by  the  medium  as  a  whole  or  by  the  walls  containing  the  medium. 

Since  explicit  inclusion  of  the  effects  of  the  walls  in  a  calculation  is  not 
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feasible,  we  will  treat  the  excess  momentum  as  belonging  to  the  medium 
as  a  whole.  From  this  point  of  view,  the  propagation  of  the  field  should  be 
considered  relative  to  the  center  of  mass  of  the  medium.  Therefore,  in¬ 
stead  of 

lo1*1  La/i exp Wr*]  +  m\  exp<-l#3,J 

(where  ag  and  s^  are  respectively  destruction  and  creation  operators, 
apart  from  constants),  we  use 

$<*>  ■  # 

=  (  Mtr r)irA.«p(fc2- j;A>  +  VAexp(- r J 

When  E(R)  is  applied  to  a  state,  it  not  only  creates  (destroys)  a  photon, 
but  also  alters  the  momenta  of  all  the  atoms  by  a  small  and  equal  amount 
in  such  a  way  that  total  momentum  is  conserved.  From  now  on,  therefore, 
when  we  write  down  a  state  it  is  to  be  understood  that,  in  addition  to  the 
photons,  excitations,  etc.  specifically  noted,  the  state  also  includes  an 
additional  superimposed  momentum  of  the  entire  medium  (equally  divided 
between  the  N  atoms)  such  that  the  total  momentum  of  the  state  is  zero. 

With  this  understanding  all  the  states  we  are  talking  about  have  the  same 
momentum,  and  we  can  study  the  propagation  of 

<°IEvn>  =  £x  <oiEvi7x7rx>. 

There  will  be  no  need  to  mention  this  explicitly  again.  Indeed,  most  of  the 
subsequent  manipulations  could  be  gone  through  formally  without  performing 
this  trick  at  all,  the  only,  trouble  being  that  they  would  be  meaningless 
due  to  both  sides  of  the  various  equations  being  zero.  The  energy  associated 
with  this  superimposed  motion  is,  of  course,  negligible  in  the  limit  N— 

We  now  give  some  typical  Hamiltonian  matrix  elements.  All  the 
ones  we  will  need  will  differ  from  these  at  most  by  the  presence  of  an 
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excess  momentum  on  an  atom  appearing  on  both  sides  of  the  matrix 
element  and  having  no  effect  except  to  shift  the  energy  In  the  case  of  a 
diagonal  matrix  element.  The  energy  of  the  state  j  0^>  Is  taken  as 
zero.  We  obtain  the  following  matrix  elements  for  the  Hamiltonian 
operator  J/  : 

<^U,/^A)/#/Al.£(A)>  =  It  1*0  + K fo* .£>]  *  (1) 

where 

K^yc)  = -fe  t2  a- 3- *  ^21  •  <*> 

;  (3) 

/  /-///V  2  U1/2  /  8  A/2  1/2 

=  •  ^hrvo/*  i(7v^}  fi-  <y);  (4) 

<M.;*l(A)l$IBi.stW,(sZ1  -£)(*)>  «  /^~^t  (5) 


In  equations  (4)  and  (5),  A  Is  the  absolute  value  of  the  dipole  moment 
matrix  element  for  the  transition,  and  G{  }C )  Is  a  cutoff  factor  which  we  now 
pause  to  discuss. 

Our  state  /  0  y  has  each  atom  In  a  definite  momentum  state,  and 
we  are  not  Including  any  Interactions  in  which  two  atoms  can  exchange  . 
momentum  with  both  remaining  in  the  ground  state;  that  is,  we  are  neglect¬ 
ing  direct  interactions  between  atoms  in  the  ground  electronic  state.  If 
these  interactions  were  included,  their  main  effect  would  be  to  introduce  a 
correlation  in  the  positions  of  the  atoms  such  that  the  distance  between  the 
centers  of  mass  of  two  atoms  can  never  be  less  than  an  atomic  diameter. 

We  take  this  into  account  by  altering  the  interaction  so  that,  while  atoms 
are  formally  permitted  to  approach  each  other  arbitrarily  closely,  they 
are  not  permitted  to  exchange  electronic  excitation  except  at  distances 
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greater  than  the  atomic  diameter  Tq.  In  this  way,  the  transfer  of  excitation 
should  proceed  just  as  if  there  were  no  pairs  of  atoms  closer  to  one  another 
than  r^.  In  momentum  space,  this  means  that  the  momentum  transfer 
accompanying  the  exchange  of  electronic  excitation  should  not  be  greater 
in  order  of  magnitude  than^l/r^),  and  this  is  embodied  in  the  factor 
G(/t  )•  Thus,  G  is  a  function  which  is  nearly  unity  when  /c  <<,  1/tq,  and 
becomes  very  small  or  zero  when  tt  For  the  Gaussian  cutoff 

factor  used  later  on  in  this  article,  a  way  of  relating  the  momentum  space 
cutoff  to  a  particular  value  of  rQ  is  given  in  the  appendix.  This  prescription 
is  not  unique,  however. 

We  will  use  the  remainder  of  this  section  to  dispose  of  one  more 
preliminary  matter.  Later  on  we  will  have  occasion  to  evaluate  the  quantity 


W  =  lim 


<  -9  0+ 


(6) 


In  situations  of  interest  to  us  here, 


Thus,  we  can  estimate  W  by  replacing  by  vq  .  With  this  replacement, 
however,  Re  W  is  just  the  second-order  transverse  self-energy  of  the  state 
I  A  a j  X  (A)^and  may  therefore  be  absorbed  into  v  by  renormalization.  If 
we  also  assume 


vo 

and  remember  that  G(vQ/c)  =  1,  then  we  get  (replacing  the  sum  over  by 
integration) 

3 

V  =  j  /**•  2  ,  (7 ) 

c 
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where  c7~  is  the  natural  linewidth.  In  what  follows,  we  will  always  use  (7). 

The  approximations  involved  are:  (1)  Neglect  of  recoil  corrections  to  the 

natural  linewidth,  which  are  small  and  can  be  included  if  desired;  (2) 

Neglect  of  of  the  dependence  of  Im  W  on  v^,  which  is  also  very  small  in 

the  frequency  region  of  interest;  (3)  Neglect  of  the  cutoff-dependent  Re  W 

when  vq.  The  formal  justification  of  this  requires  the  Arnous -Bleuler 

4  9  10 

transformation.  '  '  If  the  reader  will  accept  this  one  result,  he  will  be 
spared  the  complications  that  this  transformation  would  introduce  In  the 
subsequent  manipulations. 

III.  BASIC  ASSUMPTIONS 

Classically,  one  may  expand  the  electric  field  and  Polarization 

P  in  Fourier  series  as  follows: 
aa 

E(r,t)  =  J  dv  ^  (v)^  exp^t^-  r  -  vt]  ; 

—  CO 

P(r,t)  =  J  dv£^  F ^  (v)£c  exp  ~  *  V<1‘ 

_  -oO 

Here  the  notation  eC  represents  a  wave  number  St and  a  polarization  unit 
vector  £  s  ;  for  each  wave  number  there  are  three  polarizations,  two 
transverse  and  one  longitudinal.  The  allowed  wave  numbers  are  those 
compatible  with  the  periodic  boundary  conditions. 

The  requirements  that  the  medium  be  isotropic  and  describahle  by 
a  refractive  index  (but  with  the  possibility  of  "spatial  dispersion")  may  be 
expressed  by 

4*  Pc  M  *  •  4.C  •/*>  E.C  <v>-  <*> 

where  F  is  a  function  depending  only  on  the  frequency  v,  the  absolute  value 
/t  of  the  wave  number,  and  whether  the  mode  «£  is  longitudinal  or  trans- 

We  will  also  use  the  notations  Ft(v,  *  ),  F^v.x  )  to  indicate  the  trans 


verse. 
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verse  and  longitudinal  susceptibilities  respectively. 

Quantum  mechanically,  E  and  P  become  operators,  but  other¬ 
wise  the  situation  is  unchanged.  In  the  Schrodinger  Picture,  the  frequencies 
are  to  be  replaced  by  energy  differences  divided  by  ii,  and  we  have  different 
equations  for  the  various  matrix  elements  of  the  operators.  A  quantum 
analog  of  eq.  (8)  is 

<0  I  4.  P,-  I  X  >  =  F(V  ^4'^)  <0  |  Ex  (  X>.  (9) 

It  will  be  convenient  to  rearrange  (9)  slightly  before  using  it.  We 

note  that 

E  =  E  -  —  A 

,v\  mcouI  c  at 


where  Ecou^  i®  the  electrostatic  longitudinal  Coulomb  field,  and  A  is  the 
transverse,  vector  potential  (we  use  Coulomb  Gauge  throughout).  We  also 
have 

<°'  -4  <°U  >x>-  <io» 


Hence  if  we  define 


£  - 


E  ,  +  lV\  A, 
vi  coul  A* 


then  we  can  combine  (9),  (10),  and  (11)  to  obtain 


<ol  41.  Pj.  I  X>  «F(V  Aj.  ^.|<0  |  (12) 

Eq.  (12)  is  the  expression  we  will  use  for  the  requirement  that  the  medium 
be  describable  by  a  refractive  index. 

It  is  also  useful  to  notice  that  the  electrostatic  dipole -dipole  inter¬ 
action,  whose  matrix  element  is  given  by  (5)  may  be  written  as 

A  m*cou! 


(B,A), 
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where \  electric  dipole  moment  operator  of  atom  A  and  Ecou^ 

(B,  A)  is  the  electrostatic  field  produced  by  atom  B  at  the  position  of  atom 
A.  We  can  therefore  write 


(B),  Ig  )(A|  > 

■  <*>>■  031 

We  also  note  that  eq.  (4)  can  be  rearranged  as  follows: 

=  j 

X 

■"  -ji  i  <14> 


since  we  are  interested  only  in  the  situation  where 

I  vx  -  vQ|  <<  vo  .  (15) 

In  order  to  obtain  a  closed  set  of  integral  equations  to  determine 
the  functions  F,  we  will  also  need  a  stronger  assumption  which  Is  a  gen¬ 
eralization  of  (12): 

(16) 

<<  <A>/4*P-(A,I  >>*  -  k<sa:*>. 

(A) 

where  Pv  '  is  the  polarization  of  the  medium  due  to  all  atoms  except  atom 

yv\ 

A.  To  understand  the  meaning  of  (16),  suppose  that  the  initial  state  of  atom 
A  is  not  known  with  certainty,  so  that  the  "initial"  state  of  the  medium  is 
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not  simply  j  0  ^>  ,  but 

|0' >  =  |  0>  +  ul/fc  (A)  >  . 

Now  the  matrix  element  0'j^.Jx.^  will  have  a  term  of  frequency  - 

k<Sa>  &  )  equal  to  |^^and  a  similar  situation  will  hold  for 

(A)  **  ** 

Px  .  Equation  (16)  expresses  the  assumption  that  the  polarization  of  the 

/vv. 

rest  of  the  medium  is  still  related  to  the  field  (which  now  includes  contri¬ 
butions  from  the  virtual  excitations  of  atom  A)  by  the  same  susceptibility 
function.  Eq.  (16)  does  represent  an  approximation  beyond  (12);  Classically, 
we  can  always  limit  our  considerations  to  a  linear  susceptibility  simply  by 
making  our  fields  arbitrarily  small;  but  in  quantum  mechanics  we  cannot 
make  the  matrix  elements  appearing  in  (12)  become  small.  Their  values 
are  fixed  by  the  solution  of  the  eigenvalue  problem.  Nevertheless,  we  can 
define  a  function  F  by  means  of  (12),  and  this  will  lead  to  no  contradictions 
as  long  as  we  restrict  ouselves  to  those  matrix  elements  in  (12).  Moreover, 
this  is  the  obvious  quantum  mechanical  analog  to  the  classical  procedure  of 
allowing  the  fields  to  approach  zero:  the  one -photon  state  is  as  close  to 
zero  as  we  can  get  and  still  have  any  field  at  all.  Strict  linearity  would  be 
expressed  by  requiring  that  expressions  such  as  (12)  hold  for  all  matrix 
elements  of  the  operators  involved,  and  in  particular  for  those  appearing  in 
(16).  Hence,  (16)  is  a  stronger  linearity  assumption  than  (12).  In  this 
article  we  will  assume  that  it  holds,  but  in  the  next  article  it  will  be  shown 
how  the  failure  of  (16)  to  hold  exactly  can  lead  to  changes  in  the  value  of  the 
function  F  appearing  in  (12). 
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IV.  INTEGRAL  ECUATIONS  FOR  F 

In  order  to  find  expressions  for  the  various  matrix  elements 
appearing  in  eqs.  (12),  (16),  we  must  know  something  about  the  dressed 
photon  state  }  X.^>  .  Denoting  general  states  by  a,  b,  etc.,  we  have 

jfcl»X-\>  <•!*>  ■  4>*a<’a(#|l>><b  |t>.  (17) 

where4i  is  the  energy  of  the  zero-order  state  j  a.  ^  .  The  appropriate 
"retarded"  solution  of  (17)  is 

<a  |\>  =  ^_1(vx  -  va*  ij5>)'1^b^a<a!X|b><b  |  X>  (18) 

in  which  the  limit  ^  — >  0+  is  always  understood. 

Now  let  us  examine  some  field  and  polarization  matrix  elements. 

If  /jf  is  a  transverse  mode,  so  that  there  are  photons  of  type  ,  we 
have 


<°  l£  i  I  *>  ■  <»  I  *•  > 


(19) 


For  the  polarization  we  have 
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so  that 

<0  /  4*P^/\>  =  £asOI  4rF,y  l  Aj./^(A)><Aj,/^(A)  /  X> 

=  <Aj.^(A)  1 1>  .  (20) 

Here  and  henceforth  a  summation  over  repeated  coordinate  indices  j,  k, 
etc.  is  understood.  We  will  now  apply  (17)  to  /  a  ^  =  /  Aj,  M (A)  ^  , 
and  also  use  (14).  The  result  is 


<"j.^(A)l  x>  =  yA*p<0 /£^/yX^/x> 

+  <?st?<Ai-zr,A>IW(zr -gsH*)xs.m-tsHA)  i\> 

Ct)  a 

(®l  I ~x> 


Here  the  last  summation  goes  over  all  longitudinal  modes  «c.  •  and  we  have  ! 
defined 


^  S-  VX  -  vo 

We  now  apply  (18)  to  the  state  /  a^  ~K*>  HA)  V  ,  and  also 

AA  £  m!5 

use  eqs.  (6).  (7),  and  (19).  The  result  is 

(21) 

^[w  -/<(« lA*>5ry,]<Aj.^tA)|X>  =  X)-ior(Aj,/fc^(A)  ) 

*  <aM /&•  <*7--tc  ><*>.&  <B>XP>t- ytr4t)(A>’4c  ‘B»  I  'K>. 

The  last  summation  now  goes  over  all  modes,  transverse  and  longitudinal. 

In  arriving  at  (21),  we  have  assumed  that  the  omission  of  the  one  photon 
mode  At  from  the  summation  in  (6)  does  not  affect  the  result.  The  matrix 
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elements  of  TU  are  defined  by 

<Aj  (A  )(K  IS v  ^(B).  ( £  )(A)  >  =  ( Aj,*(A)  j#/^,  AL  (B),  {/t  ~^)(A)  > 

<(Aj ,/i(A)I^U,{k  -&)(*)><£>* -'UWWh'Jk  <BM£-£)(A)> 

+  —  (22) 

' *  -^(2a.^-4c))  +  l3 1 

The  first  term  on  the  right-hand  side  of  (22)  is  nonzero  only  for  longitudinal 
modes,  the  second  only  for  transverse  modes. 

We  now  make  the  following  definition: 

t-:r&  * 

x<Bk,l/rA.  >(*>.£  (B)|\> 

S  , /C^y  ,^^)  ^A\k.,JCj{A)  J  .  (23) 

Now  eq.  (23)  can  be  inserted  into  (21),  and  the  result  solved  for^ftj,/c^  (A)jx^ 
The  result  is 

^Aj,^(A)|x>=  -•<£  <o|£^)x>[0>  -K(%A' fy) -*]■££  2*  (24) 

The  matrix  U  -  -  K  -  A  )  *  is,  of  course,  the  inverse  of  the  matrix 

whose  elements  are 

(UJ-  K.)Sj&  -  b  i*‘ 

It  is  now  a  simple  matter  to  substitute  (24)  into  (20),  sum  over  A  using 
the  Boltzmann  distribution,  and  solve  eq.  (12)  for  with  the  aid  of  the 
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result.  This  gives  the  solution 


« 


where 


/7  g- 

O'  =  — E -  ,  and  k  is  Boltzmann's  constant. 

2  mkT 


We  could  have  replaced  the  matrix  product  under  the  integral  sign  by 


(p  Ujk£7k  = 


It  will  also  be  necessary  to  derive  an  expression  analogous  to  (25) 
for  the  longitudinal  field.  If  is  a  longitude  mode. 


<o/<£Ji>  *£a<o /<£x  I A  j,  /Cj.  (A)  X  Aj»xraC(A)  <26> 

and 

<0  /i.pjk  >.  -Sit  B  ^<6k.  ^(B)  /  i  >  (27) 

Now  eqs.  (17),  (18)  may  be  used  much  as  before  to  give 

-U  cr  (B)/ 
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By  using  (6),  (7),  (13),  (22),  (23)  in  a  manner  similar  to  the  derivation  of 
(21),  we  easily  find 


(28) 


»'[*  (B)  |X>  = 

+  *  'Sb1^  • 

in  which  we  have  also  made  the  assumption  that  the  omission  of  one  atom 
from  the  sum  in  (26)  will  not  noticeably  affect  the  result.  Now  one  proceeds 
just  as  in  the  derivation  of  (25):  Solve  (28)  for  ^B^»  (B)  I  X  insert 

the  result  in  (27),  carry  out  the  summation  (replacing  it  by  integration 

over  the  Boltzmann  distribution),  and  solve  eq.  (12)  for  F^.  The  result  is 
.2  M  A  3/2/ 


F^vJ,  fc)  -  - 


3/2 /• 

■“  < -^-)  (  exp(  ^q2 )  [t^  K(q.  ,*  )  -  A  (u>,  ^q)]^1 


(29) 


Eqs.  (25)  and  (29)  determine  the  two  functions  F  if  A  can  be 
determined.  We  must  now  turn  our  attention  to  the  evaluation  of  A  .  There 
are  (as  far  as  we  know)  two  methods  for  doing  this.  One  uses  a  formal 
infinite  series  expansion  and  is  perfectly  general  if  the  series  converges. 
The  other  gives  the  result  in  closed  form,  but  requires  the  assumption 
(16)  which  has  not  been  used  up  to  now. 

(1)  Formal  Expansion 

By  comparing  eqs.  (20a),  (21),  (23)  it  is  easy  to  see  that 


*  ’Sa*  <^Ak*  (A>  I 

=  <^a<Ak,  I  a  ><a  /  X  ^  . 


(30) 


The  summation  goes  over  all  states  except  I  a^>  =  Since  there  is  no 

singular  contribution  from  Z'*^,  it  may  be  considered  to  go  over  all  states. 

The  next  step  is  to  substitute  (18)  into  the  right-hand  side  of  (30). 
Then  use  (18)  again  to  replace  all  the  terms  ^  b  I  X)on  the  right-hand 
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si.de  of  the  result,  except  for  /  b.^  =  j  A^>  /£,  (A)  ^  ,  and  repeat 

this  indefinitely.  This  gives 

*  V^’MK  V#(A)/r> 


<Cftj,^(A)/^/a><a/^/bXbfx> 
1  b  Mvx-*a+iS*) 


(A)/3y/aX^/Ak./fc(A)><rAkt^(a)i\> 

*  (v\  -  Va  +  1  S’  > 


+  <Aj.  ^(A)/j^/aXk/j^lbXbl\> 

b’  *  K  -  va  +  1  5  ) 


a  b 


-  /^Ajt 5 (A>/^><a/j/  ^k./g  ) XAk»  ^  (A)  \  ~X.> 

H  (vx  -  va  +  i  ^  ) 

f'  <A j .  *  (A ) ijl |aX a|^|b > <bjj/ |  A  fc.  £  (A)><Ak,  g(A)  1 k> 

+  a,b  - 2 - Z> - 3 -  +-( 

*  (V  -  V  +  if  )  (v  -  V  +  i  f  ) 


where  a  prime  on  a  summation  means  that  the  three  states  /  A^,  /X  (A)^ 
are  to  be  ommitted  (k  =  x,  y,  z).  By  inspection  of  this  series,  we  see  that 
can  be  expressed  formally  as 


i.L  k>#<A» 

jk  ’*  ’”A  a  K  (vx  -  va  +  if  ) 


<Aj,  £(A)/ifa  ><  aj^/bXb/^/A.,  £  (A)  > 

& *  !  »  ■—  —  . . .  i  ■■  H  -  mi  —  -■  ■■  -  ii  ■ 

a,  b  n  -  v  +  i^Hv  -  v.  +  if  ) 


+  -—  02) 
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The  expansion  (32)  will  be  made  use  of  In  the  following  paper.  For 
the  remainder  of  this  article,  however,  we  will  make  use  of  a  formal  closed 
form  solution  which  depends  on  the  assumption  (16). 

(2)  Closed  Integral  Equation  Using  (16) 


We  first  consider  the  contribution  of  the  longitudinal  coulomb 

part  of  to  A ...  Denoting  this  by  a  superscript  H  ,  and  using  eqs.  (5), 
jk  -  * 

(22),  and  (23),  we  can  express  this  term  as 


*  Ajk(w'£  '5a) <Ak*^  <A> I  *>  - 

-  A.  wMB*- &*.  <B>-  (B).  <* -&><*>  I  \> 


where 

«(A) 

‘  ^-BtA  ~  <6  (B)’  <5  •A)(A»  rX>  ' 

and  the  sum  over  -C,  of  course,  includes  only  longitudinal  modes.  Using 
(5),  (13),  and  (33),  we  can  now  obtain 


=  <{Z'fic){A)l^ J  Ato£(A)><Ak*£(A)  J 

<bu*-£.><a)I1>. 

=  -^Gu*)  <Ak.*c(A)/l>  -  ‘ 
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It  is  also  evident  that 


<(A  )(A)  /  4.P^(A|  I  i  (35) 

pV(A\ 

It  is  now  a  simple  matter  to  insert  (34)  and  (35)  into  (16)  and  solve  for  \  ; 


j'1  £(**)FlIu  -  X  ($A’  ]  (A)  1  ^ 

^  1  +£(>£) F/u)-K(qA.f  *4).^] 

The  next  step  is  to  substitute  (36)  into  (33),  replace  the  sum  over  •<.  by 


integration,  and  extract  the  value  of  A..  •  The  result  is 


Mi.  y.  3 

■II  J _ *  A 


If  A.^(w,/6 ,  q)  -  ■  /  .  . — 

2w2  J  1  +&(y£')F,U>-J<  (q,,fc -/£'), /fc'l 

V  1  Art  *\  ™ 


There  is  one  more  contribution  to  A..  »  namely  that  due  to  the 

JK 


transverse  part  of  .  It  can  be  written 


v02/^Z6c  k  (B) 

^‘7  v  v^[vx-v^  -/f  (qA.>f-^)  +  ijl 

t  (A) 

_  _ _r  (38) 

=  v7  [vx-v7  -K(qA.^-^l  +  ‘31  ^ 


where 


£vk  <6k.(«-A^MA).A«(B)ri> 
BiA  7  M  ^7  *  * 
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and  we  have  used  (22),  (23),  plus  an  obvious  generalization  of  (4).  We  now 
proceed  much  as  in  the  derivation  of  (37).  First,  by  using  (4),  (14),  and 
(18),  we  find 


=  <o  XtfM  1  >-> 

pt  (A) 

^/<0/£^/^>/2[^^k<Ak,^(A)  I  X>  +  l7  ] 

*K*v?  + 1  5  i 


(39) 


We  also  have,  obviously. 


t(A) 


(40) 


i— 1 1  (A) 

Now  one  substitutes  (39)  and  (40)  into  (16),  solves  for  ,  inserts  the 

result  into  (38),  and  again  extracts  the  value  of  Afc  .  The  result  is 


Jl  Aj.  (K/£  .q) 

JK  MM 

4.  2 


(41) 


voA  C[  £  jk'  ~^7T~  ]  6-2(^,)Ft[w-/f^^.^«),/!g»]d3>ig» 

8irZ  J  ^  c2/^,2(vx  -  c/t'  -  /<(q,  *-/£')  +  i5)2 

+  VQ2  fi(/L')c  /t*Ft[w  -  /f  (q,*-  -/*,').  /‘t-'] 


( 


X  (v, 


-  c/t’  -  K(q,  +  iS  ) 


where  we  have  also  used  (4)  and  (14)  to  get  the  matrix  element  of  £f  . 
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The  matrix  elements  of  A  are  now  given  by 

v  *  -  ‘  ^  ^  jk +  “U  +  (42) 

Equations  (25),  (29),  (37),  (41),  and  (42)  now  determine  both  functions  F. 

The  observed  refractive  index  ^  is  determined  by  the  requirements  that, 
first,  the  electromagnetic  waves  propagating  in  the  absence  of  macroscopic 
charges  are  purely  transverse;  and  second,  that  the  usual  relation  between 
wave  length  and  susceptibility  holds.  Thus  we  write 

(  M  -  1  =  Ft(W,fVx/c)  =  Ft(W,  C^c).  (43) 

The  solution  of  (43)  will  in  general  require  the  analytic  extension  of  Ffc  into 
the  upper  half  of  the  complex  A.  -plane. 

V.  STATIC  LIMIT  AND  LOWEST  ORDER  CORRECTION 

In  the  previous  section,  a  set  of  coupled  integral  equations  was 
derived  which  in  principle  determines  the  refractive  index  if  the  assumption 
(16)  is  valid.  However,  due  to  their  highly  nonlinear  character  (among  other 
things),  it  is  evident  that  a  closed  form  solution  of  the  equations  is  not 
feasible.  Therefore,  suitable  approximations  must  be  found.  The  reader 
will  remember  that  we  are  particularly  interested  in  the  experimental  con¬ 
ditions  prevailing  in  Tomiser's  experiments  on  sodium  vapor.  It  will  be 
recalled  that  the  observed  linewidth  in  these  experiments  was  much  greater 
than  the  Doppler  width,  suggesting  that  the  effects  of  translational  motion 
of  the  absorbers  is  not  very  important;  this  corresponds  to  the  limit 
m  — >o fit  which  we  call  the  "static  limit.  "  This  limit  is  also  of  interest 
as  a  check  on  the  work  of  references  4,  5,  and  7,  in  which  the  absorbers 
were  considered  fixed  in  given  positions  in  space. 
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In  the  static  limit,  all  the  quantities  K  become  zero,  and  it  is 
easy  to  see  that  A  becomes  a  scalar: 


Because  A  is  a  scalar,  we  also  have 


x(c/fc)  x(v^  -  c/C  +  i  &  ) 

(47) 

Equations  (44)  -  (47)  are  essentially  the  same  as  the  corresponding 
equations  in  references  5  and  7.  There  are  two  main  differences:  first, 
the  use  of  the  Arnous-Bleuler  transformation  brings  about  some  small 
alterations  in  the  denominator  of  (47),  which  will  be  discussed  in  more 
detail  in  the  next  section.  Second,  in  the  earlier  work  the  cutoff  function 
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G  was  left  out  until  the  very  end  and  then  inserted  only  when  it  was  needed 
to  prevent  a  divergence;  hence,  G  does  not  appear  in  the  analog  of  (47)  at 
all,  nor  in  the  denominator  of  the  analog  of  (46).  Neither  of  these  differences 
is  very  significant;  the  omission  of  the  Arnous -Bleuler  transformation  and 
the  way  of  bringing  in  the  cutoff  are  both  rather  arbitrary,  nonessential 
elements  in  the  formalism.  Since  Fg  is  a  function  only  of  u>,  we  have  in 
this  limit  CX(UJ)- 1  =  Fg  (**>). 

We  would  now  like  to  discuss  the  lowest  order  finite-mass  cor¬ 
rection  to  the  static  approximation,  but  this  will  be  made  easier  if  we  first 
make  another  approximation.  A  comparison  of  (37)  with  (41)  and  (7)  shows 
that 

>  >  <s~.  . 


This  is  because  only  the  integral  (37)  requires  the  cutoff  for  convergence, 
making  t£  proportional  to  the  third  power  of  the  cutoff  wave  number. 
For  the  remainder  of  this  section,  therefore,  we  will  set  A  equal  to  it  ; 
we  will  also  omit  the  G  in  the  denominator  of  (37)  since,  as  pointed  out 
above,  the  way  in  which  the  cutoff  is  brought  in  is  rather  arbitrary. 
Equations  (29),  (37)  now  become 

(48) 

F1(W,/C)=  -  D(P/*)  f  exp(-Pq2)^[u)-/<(q,^)-A(U,^]^  ^  d3q; 


A.. 

JR  *  % 


(49) 


where 


4ir>i2 

K 


B 


•■4- 

2i rn 


(50) 


9 
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We  now  expand  and  in  inverse  poweres  of  m,  with  the 

intention  of  evaluating  the  zero-  and  first-order  terms.  We  write 


F/W./t,)  =  F  (U)  +  +  — 

*  8 

m 


and 


V iu’X-$  =  “<“>4  *  g  'll-  ♦  - 

'  '  m 


(SI) 


(52) 


It  will  now  be  convenient  to  use  a  specific  form  for  the  cutoff  function.  The 
choice  made  in  this  section  is 

G(  /£.  )  =  exp{-^2  yj") 

According  to  the  prescription  given  in  the  appendix,  this  corresponds  to  a 

r  1/2 

cutoff  radius  of  r^  s  * 

We  now  insert  (51),  (52)  into  (48)  and  (49),  expand  in  powers  of 
(1/m),  and  equate  coefficients  of  like  powers.  By  equating  the  zero-order 
coefficients  we  get 


F8(C)  =  - 


U-  Ag(w) 


and 


= 


-AF  s^ 


1  +  Fs(c*>) 


where 


/U  -5L 


3/2 


2^C 


(53) 


(54) 


Equations  (53),  (54)  can  easily  be  solved,  giving  the  result 


Fg(W) 


-(CJ4  D)  +  [(w  -  P)2  -  4D  A 1 

2(U>  +  A  ) 


1/2 


(55) 
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According  to  (55),  Fg  has  an  imaginary  (absorptive)  part  in  the  region 

}U  -  D  j  <  2  (D  A  )1/2 


In  this  region,  of  course,  the  square  root  is  taken  so  that  the  imaginary 
part  is  positive.  In  the  regions  to  the  left  and  right  of  the  absorptive 
region,  it  is  seen  by  analytic  continuation,  or  by  physical  requirements, 
that  one  should  take  the  negative  square  root  on  the  left,  the  positive  one 
on  the  right. 

Since  is  of  the  order  of  an  atomic  radius,  in  a  dilute  gas 


A  »  D  • 


Therefore,  the  denominator  in  (55)  is  a  slowly  varying  function  of  in  the 
absorptive  region.  If  the  denominator  is  considered  approximately  constant, 


we  obtain  for  the  half -width  at  half-height 


*  (3D/t)l/2  =  (^)1/4^£-  (-£,>* 


Thus,  the  half -width  should  be  proportional  approximately  to  the  square 
root  of  the  density  of  the  gas. 


'When  we  equate  coefficients  of  the  first  power  of  (1/m)  in  (49),  we 

get 


(57) 
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where 


U')  g  . 


(58) 


la  arriving  at  (  57) ,  we  have  made  use  of  eq.  (2)  for  K.  The  prime  on 

F  denotes  diff  erenti at  i  on  wi  th  respect  to  .  We  can  now  insert  (! 
8 

and  (57)  into  (48),  expand,  and  equate  coefficients  of  (1/m): 


g|W.  A?)  =  a(V)  +  b (O)*2  +  c(u>)  (<d). 


where 

a(O)  =  - 

b(w)  =  - 

c(W)  = 


D  B  Fg  (L>) 


3/2 


[W  -  A  (<*»)]  [1  +  F  (U)] 

a  a 


,12  8X‘ 


'it* 


D  * 


r  i 


B  F  (U) 


(-^-7) 


3/2 


2[d;  -  As(W)]2  £  3[1  +  Fa(W)]2  2<^ 


D  B 


3[  i  +  Fs(w)r 


To  evaluate  ^(O) ,  we  substitute  (59)  into  (58),  getting 
3  (W)  ^GrZ(/V)[ a(w)+  b(^)>fc'2  +  C(W)  5(^)1  d3/t' 

*  (  —  ,)3/2  [a(W)  +  +  c(^)  5  («)  ]  . 

2-C 


(59) 


(60) 


(61) 
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Equation  (61)  i  s  now  solved  for  J*  (  uJ  ): 

r  ,  *  ,3/W>  ♦  -Lm-  ] 

O  M  =  2oC C _ 4^. 

,  ,  it  v  3/2 

1  _  (  — -  .  )  c(<4) 

2£T 


(62) 


Equations  (59)>  (60),  and  (62)  determine  g.  Since  A  is  still  a  scalar 
in  this  approximation,  we  still  have  Ffc  =  F  ^  .  Note  that  the  result  is  still 
independent  of  T. 

We  have  calculated  Im  Ft(  CJ  •  vq/c)  both  zero-  and  first-order 
under  the  conditions  of  one  of  Tomiser's  experiments  on  sodium.  F^  is 
sufficiently  small  in  this  case  that  we  can  neglect  its  dependence  on  wave 
number.  The  following  parameters  were  used: 

T  =  688.1°  K. 

P  =  0.  503  mm  Hg. 

«C  =  2.  83  x  10  ®  cm 
JjJ"  -  4.167  x  10  ^  erg  cm"* 

O'  =  3.199  x  107  radians/ sec 

\  =  2tt/v  =  5.89  x  10  ^  cm. 
o 

The  results  are  shown  in  figure  1.  Obviously,  the  lowest-order  correction 
for  finite  mass  has  only  a  very  small  effect.  If  one  looks  closely  enough 
at  the  figure,  one  can  see  that  the  finite -mass  correction  slightly  raises 
the  maximum  of  the  line;  it  is  therefore  to  be  thought  of  as  a  narrowing 
effect.  This  can  be  understood  qualitatively  as  follows:  In  the  static  approx¬ 
imation,  the  line  is  broadened  due  to  the  transfer  of  excitation  from  one 
absorber  to  another.  This  transfer  may  be  accompanied  by  a  momentum 
transfer  of  any  size  up  to  the  cutoff  momentum.  When  the  finite  mass  is 
taken  into  account,  however,  excitation  transfer  accompanied  by  large 
momentum  transfer  is  inhibited  due  to  the  recoil  of  the  atoms  spoiling  the 
energy  balance;  i.e. ,  excitation  transfer  with  large  momentum  transfer  is 
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no  longer  a  "resonance"  process.  This  effect  might  be  called  "recoil 
narrow  lag". 

The  half  width  at  half  maximum  from  fig.  1  Is  (In  wave  length 
units)  1.48  A.  The  experimental  values  are  1.  565  A  for  the  3/2  line,  and 
1.  075  for  the  1/2.  Thus,  we  are  able  to  get  qualitative  agreement  with 
experiment  with  a  reasonable  value  for  oC  .  In  the  next  section,  we  take  Into 
account  explicitly  the  doublet  character  of  the  line,  making  possible  a  more 
detailed  comparison  with  the  experimental  results. 

VI.  THE  DOUBLET  CASE  IN  STATIC  APPROXIMATION 

The  sodium  D-llne  Is  not  actually  a  transition  from  a  nondegenerate 

s -level  to  a  triply  degenerate  p-level,  as  assumed  up  to  now.  The  transition 

2  2  2 
Is  actually  from  S  to  P,  and  the  excited  state  has  two  components, 

and  ^P3/2*  The  wave  lengths  of  the  transitions  are  respectively  5896  A 
o 

and  5890  A  . 

A  detailed  derivation  of  the  formulas  for  the  doublet  case  would  be 
largely  repltltlous.  We  will  therefore  content  ourselves  with  rearranging 
some  of  the  formulas  of  the  previous  section  for  the  static  limit,  and  then 
giving  without  proof  the  corresponding  generalization  to  the  doublet  case. 

The  justification  for  the  use  of  the  static  limit  Is  the  smallness  of  the  lowest 
order  correction  found  In  the  previous  section. 

In  the  static  limit  we  have  (cf.  (53)) 

f  2  -  1  =  F  =  -  - ^ -  (63) 

V  3  U  -  A 

(We  can  drop  the  subscript  s  from  now  on,  since  In  this  section  we  are 
always  dealing  with  the  static  limit.  ) 

If  we  omit  the  G  in  the  denominator  of  (46),  andi4*e  (50)  and  (63), 
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we  find 


A*  = 


b, i - r  ) 


B  1  -  f2 

3  T“ 


)  <  — =->) 

2£- 


3/2 


2  3 

(St)  d/L 


(64) 


2  2 

where  again  G  =  exp (~JZ  /L  )•  If  instead  we  use  a  sharp  cutoff  in  config¬ 


uration  space  at  a  distance  r«,  we  find 


i 

A  = 


4ir2  B 


3  /t 


( 


fS 


) 


(64a) 


f2 

In  case  of  (47),  we  will  have  to  rearrange  the  denominator  some¬ 
what  to  get  the  same  result  which  was  obtained  in  references  5  and  7  with  the 
aid  of  the  .Arnous-Bleuler  transformation.  This  appears  here  to  entail 
approximations,  but  it  is  really  more  accurate  than  direct  evaluation  of  (47) 
would  be.  First,  we  omit  all  factors  of  G  from  (47),  since  they  are  not 
needed  to  ensure  convergence.  With  the  notation  =  v,  we  can  write  the 
denominator  as 


v(vx  -  « 


+  1$)  )  +  -y-  (  f*2  -  1) 


Next  we  replace  by  vq  ;  it  is  further  assumed  that  the  main  contribution 
to  the  integral  is  from  the  region  v  vq,  so  that  v  can  be  replaced  by 
l/2(v  +  vQ).  The  denominator  now  becomes 

-j[vZ  -  (vq  +  iS)2]  [v2  -  (°2vo2]  . 

With  all  these  rearrangements,  the  integral  is  easily  evaluated: 


•l<r  (f  -  1) 


(65) 
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Hence,  by  combining  (45),  (64a),  and  (65),  we  obtain 


A  *  I  f°-  *  -5^4  <  1  ~a  f 2-  > 

l  3 v  f 2 

Equations  (63),  (66)  can  be  solved  for  ^  for  any  value  of  w  . 

In  the  doublet  case  the  procedure  is  similar,  though  it  is  complicated 

by  the  presence  of  the  spin-orbit  splitting.  One  may  use  either  a  simple 

generalization  of  the  methods  of  references  4  and  5,  or  of  reference  7,  or  of 

the  earlier  sections  of  this  article.  (Actually,  we  used  the  methods  of 

2  2 

references  4  and  5).  We  denote  the  frequencies  of  the  Pjy2  an<*  ^3/2 
transitions  by  vj  and  v^>  respectively.  Instead  of  (63),  we  find 


f2_!  .  .  °f _ i _ 

r  »K-v- 


\  -  v2  -  4 


The  equation  for  a  remains  formally  unchanged:  (66)  still  holds,  with  B 

given  by  (50). is  now  understood  as  the  orbital  dipole  matrix  element  in 

which  spin  effects  are  neglected.  We  have  neglected  v\  ~  v2  against  v^, 

so  the  natural  linewidth  is  treated  as  the  same  for  both  lines. 

We  have  calculated  the  absorption  curves  for  a  number  of  points  in 

o 

the  region  of  Tomiser's  experiments,  using  eqs.  (66)  and  (67)  with  Xq  ~  8  A, 
a  value  which  is  reasonable  and  seems  to  give  good  agreement. 

In  carrying  out  these  calculations,  it  is  convenient  to  define 

-o  ■  1/2  <"l  +  »2>  • 


r  -  ^2-v-1- , 

p  _  6jt_c^  N 

•  v  3 


Figure  2  shows  a  typical  absorption  curve,  corresponding  to 
P  =  109.  78  (T  =  688. 1°K,  P  =  }0.  503  mm  Hg. )  iXhe  comparison  between 
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experimental  and  theoretical  linewidths  is  plotted  in  figure  3.  It  will  be 
observed  that  there  is  qualitative  agreement  for  both  components  of  the 
doublet,  even  though  the  theoretical  widths  are  approximately  proportional 
to  the  square  root  of  the  density  (cf  eq.  (56)),  while  the  experimental  depen¬ 
dence  appears  to  be  linear.  The  theoretical  widths  obtained  by  the  methods 
of  reference  1  would  be  too  small  to  be  seen  on  a  graph  of  the  scale  of  fig. 

3,  so  we  are  justified  in  claiming  to  have  improved  the  agreement.  In  the 
following  article,  it  will  be  shown  that  the  discrepancy  which  remains  may 
well  be  due  to  the  failure  of  the  approximation  (16). 


APPENDIX 


We  can  introduce  a  cutoff  into  the  ordinary  electrostatic  interaction 
between  two  charge  densities  ^  \  anc^  ^  by  writing 

^2(l2)^(r12)d3r’d3r2 

s  ^2)exptiM*(5l  -*2)ld3rld3r2 

=  )exP<L^'ri)d3rij^2^2)exP<-1VS2)d3r2  j  , 


(A  -  1) 


Here  the  sum  goes  formally  over  all  satisfying  periodic  boundary  con- 
ditions  in  the  volume  V.  If  G  =  1,  we  have  the  usual  electrostatic  interaction. 
If,  both  charge  distributions  can  be  treated  as  dipoles,  we  make  the  replace- 


(A-2) 
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The  matrix  elements  corresponding  to  (A-2)  are  given  by  eq.  (5).  Now,  we 
can  determine  an  appropriate  choice  of  G  just  as  well  by  looking  at  its  effect  on 
(A-l)  as  by  examining  (A-2).  In  other  words,  we  consider  the  effect  of  G  on 
the  function 


e*p(  % 


(A-3) 


The  criterion  we  use  for  relating  G  to  a  cutoff  radius  rg  is 


lim 

r 

r  — 4  -CA 


(A -4) 


and 

4?  (0)  =  I / r Q  .  (A-5) 

Thus  i®  to  become  constant  as  r  becomes  <■<.  rg,  causing  in 

(A-2)  to  become  small;  it  is  to  behave  as  an  ordinary  electrostatic  potential 

2  2 

for  r  »rg.  If  one  inserts  G  *  exp(-,£  /c  )  into  (A-3),  one  finds  (replacing 
summation  by  integration,  of  course) 


From  (A-5)  and  (A -6),  we  find 

rQ  =  O ri/Z  ,  (A-7) 

which  is  the  relation  between  the  wave  number  and  configuration  space 
cutoffs  mentioned  in  sec.  Y.  This  prescription  will  also  be  used  in  the 
following  article.  It  is  obviously  not  unique,  however.  For  example,  if 
we  required  consistency  between  (64)  and  (64a),  we  would  get  a  different 


result. 
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